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Abstract 

We consider a nonrelativistic Chern-Simons theory of planar matter fields 
interacting with the Chern-Simons gauge field in a SU(N) g i b a i X U(l)i oca i 
invariant fashion. We find that this model admits static zero-energy self- 
dual soliton solutions. We also present a set of exact soliton solutions. 
The exact time-dependent solutions are also obtained, when this model is 

considered in the background of an external uniform magnetic field. 
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Soliton solutions in Chern-Simons ( CS ) gauge theories have received considerable 
attention over the past few years due to their possible relevance to the planar condensed 
matter systems. It is known that the abelian Higgs model with a CS term admits finite 
energy charged vortex solutions PJ • Further, the pure CS Higgs theory admits static self- 
dual soliton solutions with a </> 6 -type scalar potential 0. Moreover, in the nonrelativistic 
limit of this theory M, the charge density solves Liouville equation at the self-dual limit, 
all of whose solutions are well known. When this nonrelativistic model is modified by 
including an external magnetic field [|j] or a harmonic force ||, exact time-dependent 
soliton solutions can be obtained. The self-dual nonrelativistic case for the nonabelian 
gauge group has also been considered , which provides a unified dynamical framework 
for a variety of two-dimensional nonlinear equations M. 

In this Brief Report, we consider a nonrelativistic CS theory with a gauge group as 
in the case of semilocal Nielsen-Olesen strings |J or semilocal charged vortices |9||ll| • In- 



particular, we consider the Jackiw-Pi ( JP ) model H but with the gauge group enlarged 
to SU(N)gi b a i x U(l)i oca i. We find that this model admits static zero-energy self-dual 
soliton solutions. Interestingly enough, we are also able to find a set of exact soliton 
solutions. These solitons are characterized by the magnetic flux $ = — + 1)| n |, 

the charge Q = — -$ and the angular momentum J = Q, where n is the winding number 
and k and e are two dimensional constants to be discussed below. We also present exact 
time-dependent solutions of the model in the presence of an external uniform magnetic 
field. 

Consider the nonrelativistic Lagrangian 

£ = i*t ( dt + ieA °) V - -L| fa + ieA t )m | 2 + 9 -{^mf + ^ a A^F va (1) 
v ' 2m 2 4 

where \1> is component scalar field, i.e., ^ = (^q, -0*, . . . , V'tv-i) ( Here * denotes the 

complex conjugation ). The Lagrangian (|l|) is invariant under a SU(N)gi j, a i x U(l)i oca i 



transformations. For N = 1, the Lagrangian ([J) essentially describes the JP model. The 



N = 2 case was previously discussed and some exact solutions were obtained in Ref. ||11|| . 
Note that the scalar field self-interaction may be attractive or repulsive according as g 
is positive or negative respectively. However, as we will see shortly, the self-interaction 
is always attractive for zero-energy self-dual soliton solutions as in the case of JP model. 
The equations of motion which follow from (ffl) are 



K 



vaf) T? 



(2) 



id t m = — —DiD-9 + eA ^ - g\ V 
2m 



where the conserved matter current J u is given by, 



(3) 



(4) 



The zero component of @, i.e., the Gauss law implies that the solution with charge 
Q also carries magnetic flux $ = —-Q. The Eq. (|3|) is a 2 + 1 dimensional gauged 
nonlinear Schrddinger equation where the gauge-field variables can be expressed solely 
in terms of the matter- field variables with the help of Eq. ([|). 
The energy for the Lagrangian (|T]) is 



E 



d x 



1 



(5) 



which can be rewritten using the Bogomol'nyi |12j trick as 

1 



E 



d x 
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where a surface term has been dropped since it vanishes for the well behaved field- 

2 

variables. Now note that for the choice of g as g = =F— , the energy satisfies the bound 
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E > 0. The bound is saturated when the following first order self-dual equations are 
satisfied, 

(D 1 ±iD 2 )^ = 0. (7) 

It should be noted that the Eqs. ([7p are identical to the corresponding equation of Ref. 
except that \l/ now is a N component scalar field. 

In order to solve Eqs. ([I]), we write down the gauge potential A{ in the Coulomb 

gauge, 

Ai = -- e e l3 d jX ■ (8) 
Now it is trivial to check that Eq. (0) can be rewritten as, 

(di ± id 2 )e Tx ip j = . (9) 
Thus we have the general solution to Eq. (H) in the form ]13| , 



^ = e^f^z) , (10) 

where z = x + iy and /y(z)'s are arbitrary analytic functions. With the help of the 
Gauss law and Eqs. @ and fllCf), the decoupled equation for the x is 

P 2 N-l 

V 2 X = -e ±2 *Y,\m\ 2 (11) 
i=o 



where V 2 = 49 2 9j. Note that Eq. (|TI|) reduces to the Liouville equation in case the 
summation on the right hand side is equal to some real constant. However, we are 
interested here in more general solutions. Let us first discuss the simplest case of iV = 2, 
in which case Eq. fllTD takes the form i.e. 



V 2 X = -e ±2x (|/o| 2 + |/iW| 2 ). (12) 



The above equation can be solved exactly provided we assume a particular form for f\ (z) 
in terms of fo(z). In particular, we choose 



h{z) | 2 = | f (z) | 2 I / f (z)dz I . 



(13) 



The solution of Eq. ( |12|) is 



X = T^ln 



where sign of k in (0) must be opposite to that of ± and a 



\k\- 



(14) 



Note that for this 



choice of sign g = is always positive, and hence the scalar field self-interaction is 
attractive. If one now chooses fo{z) = ncQZ n ~ 1 {z = re ld , | n | > 1) then one obtains the 
rotationally symmetric solutions 

\/6| n I / r \ f 



e i(|nb'+|n|-l)0 • =01 



(15) 



y/ar \r , 

where 8 — \ n \(2j — 1). 

As far as we are aware off, for arbitrary iV no exact analytic solution of the equation 
(|Tl"|) is known and at present we don't know how to solve Eq. ([11]) exactly except in few 
specific cases. For example, the equation (|TT|) can be solved exactly if one assumes 



|2 JV-1 



Cj I fo(z) | 2 | / fo(z)dz 



(16) 



where N 1 Cj 



j\(N-l-j)\ 



. After substituting (|I6| ) into (p|) , we have 



V 2 X = -e ±2x |/oW I 2 
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fo{z)dz 



2\ N-l 



As in iV = 2 case we fix the convention that the sign of k must be opposite to that of 
±. Now one can check that 



N+l. 
X = =F — - — In 



a(l + \ J f (z)dz 



(18) 



solves Eq. (17), where a 



L2(AT+1)J 



- 1 . The ifij's are thus given by 



N-l 



%N+1) f (z)\Jfo(z)dz] 



a 



N+l 



(19) 



l + \Jf (z)dz 

It should be mentioned at this point that the familiar Liouville solution can be embeded 
into the SU(N) g i b a i x U(l)i oca i invariant theory for any N by choosing all the /j(z)'s 
equal. Further, for any N' < N, the solutions as given by fll8|) and (|19]) can be embeded 
into the higher N theory. 

The radially symmetric solutions for ^'s can be obtained from Eq. (|l9|) by putting 
f(z) = ncQZ n ~ 1 (z = re %e ). We find 



(.N+l) 



J(\n\j+\n\-l) 



■ , v , ■ , (20) 
v ol r V ro J L V ro J \r J J 

where 5 = \ n \(2j + 1 — N). Note that the single valuedness of ^-'s demands that 
| n | necessarily be an integer. All the p/s (pj = ipjipj) vanish at asymptotic infinity as 
r -2(|n|A r -|n|j+i)^ i m piyi n g that the rate of fall off is higher for lower values of j, reaching 
a maximum at j = 0. Near the origin p^s behave as r 2 ^ n ^ +1 ^~ 1 ^ so that all the p^s 
are nonsingular except po, which is nonsingular only when | n \ > 1. We shall therefore 
restrict ourselves to | n \ > 1, throughout this paper. 
The gauge potential A 2 is given by, 



An 



n \{N+1) 



er 



1 + 



2|/i 



-1 



(21) 



Near the origin the gauge potential goes to zero as r 2 ' 71 ' and interestingly enough it is 
independent of N. However, at the asymptotic infinity er^ ~ — | n \ (N + 1), keeping 
track of the global group structure. Infact, the profile of the gauge potential as given in 
(pl|) is same to the corresponding N = 1 case except for an overall multiplication factor 
of |(iV + 1). As a consequence the magnetic field B also has the same profile as in the 
case of JP model except for a overall multiplication factor i.e. 
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. 2n 2 (N + 1) /r \ 2 M 
B = i — ' 1 



er 2 v r 



1 + 



2M1- 2 



(22) 



so that the magnetic flux is modified to $ = — ^fy^yl n \ (N + 1). Note that the flux $ 
is N dependent and | ^$ | is quantized in terms of (N + 1)| n \. So, the flux quantum 
increases as one considers the higher values of N, i.e., enlarges the global symmetry. 
Also the flux quantum is even for odd N, while it can be both even and odd for even 
N. The charge Q and the angular momentum J are also quantized in this they 
are related to the magnetic flux by Q = J = — -$ = — | n \(N + 1). 

So far we have discussed only a set of specific solutions of Eqs. (|11|), which are 
expressed in terms of one unknown function fo(z). However, one would like to know the 
more general solutions of ([11]). Though we do not know the most general solutions of 
the Eqs. (|TT1), one can obtain a set of exact solutions in case iV = ^N'(N' — 1) ( where 
N' > 3 , i.e. N=3, 6, 10, 15, . . . ) in terms of N' unknown functions. For example, when 
N = N' = 3 the solution is given by, 



X = Tin 



where fj(zYs are chosen as 



2 



a (\M*) \ 2 + \Mz) I 2 + I0sW I 2 ) 



(23) 



fi( z ) = \tijkW jk {z) 

W i3 {z) = <pj{z)d z <f>i(z) - i,j, k = 1,2,3. (24) 



The restriction on the solution (P3|) is that the analytic functions 0j(z)'s have no common 
zeros and arbitrary otherwise. Hence in this case no rotationally symmetric solution is 
possible. However if one assumes any one of the three <pj(z) 7 s is equal to unity, then it 
is possible to have rotationally symmetric solution analogous to Eq. fll9|). Notice from 
Eqs. (PI), (|23|) and (|4]) that there is a freedom in choosing /i(-z)'s in terms of Wik(z) 
as the requirement to have exact solution is, 



2 El^-W l 2 = EU(*) 

ij=l i=l 



(25) 



The particular choice in Eq. fl2~3| ) is for notational convenience. Similar solution can also 
be written down for other values of N (6, 10, 15, ... ) using the identity, 

-2 

(26) 

where Wij(z) is defined as in Eq. (|2~3|), but now with i, j = 1,2, ... , N'. Note that the 
first sum on the right side of Eq. (26) contains N = ^N'(N' — 1) number of terms of 





" N' 




" N' 


E 1 w^z) | 2 




E 1 ^ra{z) | 2 








m=l 



the form | Wy(z) 



Let us now discuss time-dependent solutions of Eq. (|TD in case it is considered in 
the background of a uniform magnetic field. To this end notice that the action (j^) is 
invariant under dilation 



x = nr 1 *, t-^t' = n~H, $ -> = n^(t, x), 

A fc -> A' fc = 04 fc , A —> A' = Q 2 A , 



(27) 



where f2 is a constant. However, when the action (|1|) is considered in the background of 
an external magnetic field £>, only the Hamiltonian remnains a conserved quantity. This 
fact was utilized in Ref. j| to construct time-dependent solutions for JP model (N=l) 
in presence of B by starting from the static soliton solutions with B = 0. We find that 
the same conclusions are also valid for arbitrary N. In particular, the Lagrangian ([[]) in 
the presence of an external uniform magnetic field can be written as, 

C = i& (Bt + zei°) * _ J_| + ieA t - ea^ \ 2 + |(^) 2 + -^ va A^ (28) 

where Oj = —^e^Xj and = One can easily check that under the following 

transformations ( with w = — ), 
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2 wt, 
t = —tan{ — 

w V 2 ' 



1 

-tan(f 



tan(f) 



wt 

cos(—)exp 
dx v 



mw 2 /Wt, 

z r tanl — 

4 v 2 ' 



lf)j(t, x) 



i4„(t,x) 



(29) 



S = J Cd 3 x is transformed into the action without the magnetic field S = J Cd 3 x where 
L is as given by (H). Thus Eq. (29) is not a symmetry transformation of the action 
S. However, it relates the soliton solutions of ([!]) to that of the (|28|). Using the exact 
rotationally symmetric soliton solutions of the Lagrangian (fj) as given in Eqs. (^) and 
(piP, it is now straightforward to write down the time-dependent soliton solutions for 
the field- variables ipj, A and A\ with the help of Eq. (|29|) . The whole analysis goes 
through even when the solitons of (|l|) are considered in the background of a harmonic 
force 

Finally, following comments are in order: 

(i) For N = 2, the Lagrangian studied here can be obtained by taking the non- 
relativistic limit of a relativistic semi- local theory considered in Ref. ||. Infact, the 
Lagrangian of |J but with the SU(N) g i oba i x U(l)i oca i symmetry reduces to Eq. ([[]) in 
the norelativistic limit. 

(ii) It is known that the relativistic N = 2 theory admits semilocal topological as well 
nontopological soliton solutions |5| ( Actually this is also true for arbitrary iV ). However, 
as we have shown, in the nonrelativistic limit only semilocal nontopological vortices are 
admissible. Can one extend the model ([I]) and also obtain semilocal topological solitons?. 
We have checked that if the term evoA is added to the Lagrangian (0) and potential 
_|(^rt^r)2 j g m odified to |(^^ — Vq) 2 where Vq is a constant, then semilocal self-dual 



topological soliton solutions can be obtained for g = ±£- by following the discussion of 
Ref. |15|]. These self-dual solutions are characterized by nonzero energy E = unlike 
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the semilocal nontopological solitons. Further, one finds that the decoupled equations 
for the matter fields are 



2e 2 N ~ x 

V 2 lnp j = ±—(Y1 Pi-vo), j = 0,l,2,...,JV-l. 

K 1=0 

For the special case of N = 2 ( and when the constant ±— on the right side of 



(30) 



is positive ) these equations are identical to those obtained in [16,17] , and hence their 
analysis about the solutions [I7| as well as the stability (TB| goes through in this case. 



(iii) Recently, Knecht et. al. |18| have done the Painleve' analysis of JP model and 
have shown that the model is not integrable, although it naturally admits integrable 
reductions which are the familiar Liouville and 1 + 1 nonlinear Schrodinger equations. 
It would be interesting to repeat the same exercise for the SU(N) g i b a i x U{l)i oca i case. 
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